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We propose a practical folding model to describe 3 He elastic scattering. In the model, 3 He optical 
potentials are constructed by making the folding procedure twice. First the nucleon-target potential 
is evaluated by folding the Melbourne p-matrix with the target density and localizing the nonlocal 
folding potential with the Brieva-Rook method, and second the resulting local nucleon-target po¬ 
tential is folded with the 3 He density. This double single-folding model well describes 3 He elastic 
scattering from 58 Ni and 208 Pb targets in a wide incident-energy range from 30 MeV/nucleon to 150 
MeV/nucleon with no adjustable parameter. Spin-orbit force effects on differential cross sections are 
found to be appreciable only at higher incident energies such as 150 MeV/nucleon. Three-nucleon 
breakup effects of 3 He are investigated with the continuum discretized coupled-channels method 
and are found to be appreciable only at lower incident energies around 40 MeV/nucleon. Effects of 
knock-on exchange processes are also analyzed. 

PACS numbers: 25.55.Ci, 24.10.Eq, 24.10.Ht 


I. INTRODUCTION 

Microscopic understanding of nucleon nucleus (NA) 
and nucleus-nucleus (AA) optical potentials is a goal 
of nuclear physics. The optical potentials are not only 
essential quantities to describe the elastic scattering but 
also key inputs in calculations of the distorted-wave Born 
approximation (DWBA) and the continuum discretized 
coupled-channels method (CDCC) [1-3] to describe in¬ 
elastic scattering, breakup and transfer reactions. 

The p-matrix folding model is a method of deriving 
the optical potentials microscopically. Here the p-matrix 
is an effective nucleon-nucleon (NN) interaction in nu¬ 
clear matter and depends on the density p of nuclear 
matter: g = g(p). In the folding model, the optical po¬ 
tential is obtained by folding the p-matrix [4-14] with 
the target density px for NA scattering and with px and 
the projectile one pp for AA scattering; see for exam¬ 
ple Refs. [15-19] for the folding procedure. The model is 
now called the single-folding (SF) model for NA scatter¬ 
ing and the double-folding (DF) model for AA scattering. 
For AA scattering, occasionally, the optical potential is 
obtained semi-microscopically by folding the phenomeno¬ 
logical nucleon-target potential with pp. This procedure 
is also called the SF model. 

For NA elastic scattering, the SF model based on the 
Melbourne g -matrix [11] well reproduces the experimen¬ 
tal data with no adjustable parameter. In the folding 
procedure, the value of p in p(p) is assumed to be a value 
of px at the midpoint r m of interacting two nucleons: 
p = px (t* m ) • This local-density approximation seems to 
be reasonable because of the success of the Melbourne 
p-matrix folding model. 

The NA potential thus obtained is non-local, because 
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knock-on exchange processes are taken into account in 
the folding procedure. However, it can be localized by 
the Brieva-Rook method [6] with good accuracy. The 
validity of the method is shown in Refs. [20, 21]. In 
fact, the local version of the folding potential agrees with 
the phenomenological NA optical potentials determined 
from the NA scattering data [22-25] particularly in the 
surface region important for the scattering [26]. 

The multiple NN collision series in AA scattering [27] 
is more complicated than in NA scattering [28, 29]. This 
makes microscopic description of AA scattering more dif¬ 
ficult. In fact, the DF model has a basic problem. In 
nuclear matter calculations, in principle, the p-matrix 
should be derived by solving the scattering between a nu¬ 
cleon in a Fermi sphere and a nucleon in another Fermi 
sphere [30], but in practice it is evaluated by solving nu¬ 
cleon scattering on a single Fermi sphere. In addition 
to the single-sphere approximation, furthermore, in the 
folding procedure the density p of the single Fermi sphere 
is assumed to be identical with the sum of pp and px at 
the midpoint r m of interacting two nucleons: 

P = pp(rm) +PT(r ln ). ( 1 ) 

The prescription is called the frozen-density approxima¬ 
tion (FDA). The DF model with the FDA is referred to 
as the DF-FDA model in this paper. 

In actual DF-FDA calculations, the magnitude of the 
folding potential is usually adjusted to the experimen¬ 
tal data. It is then an important subject in future to 
clarify how good the DF-FDA model is. As a successful 
example, the DF-FDA model based on the Melbourne p- 
matrix well reproduced measured total reaction cross sec¬ 
tions ctr for 12 C scattering from stable nuclei at around 
250 MeV/nucleon with no adjustable parameter [19, 31- 
35]. The DF-FDA model was then applied to measured 
or [32, 33] for neutron-rich Ne and Mg isotopes [19, 31- 
35]. This analysis leads to the result that 31 Ne and 37 Mg 
are deformed halo nuclei. 
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As an alternative approach to the DF-FDA model, we 
can consider the model Hamiltonian 

H e ff = Kr + Y] Uj t + hp, (2) 

ieP 

where Kr stands for the kinetic energy with respect to 
the relative coordinate R between a projectile (P) and 
a target (T) and hp is the intrinsic Hamiltonian of P. 
Here Ui t represents the interaction between T and the 
ith nucleon in P. As Ui t, the phenomenological nucleon- 
target optical potential was used so far; see for exam¬ 
ple Refs. [1-3] for deuteron scattering. In this paper, 
meanwhile, Ufr is constructed microscopically by fold¬ 
ing the Melbourne ^-matrix with px and is localized by 
the Brieva-Rook method. The potential t/,;x is then al¬ 
ways obtainable even if no experimental data is available 
for nucleon-target scattering of interest. This is an ad¬ 
vantage of the present approach from the previous one. 
Furthermore, this approach is consistent with the fact 
that a single Fermi sphere is considered in the p-matrix 
calculation for nuclear matter. In addition, one can treat 
projectile breakup with CDCC if necessary, since Un is 
localized. These are advantages of the present approach 
from the DF-FDA model. The model Hamiltonian (2) 
can be derived from the many-body Hamiltonian with 
reasonable approximations, as shown later in Sec. II. 

It is known that the model Hamiltonian (2) well ac¬ 
counts for deuteron scattering, when deuteron breakup 
effects are properly taken into account with CDCC [1 
3]. For 4 He scattering, meanwhile, projectile-excitation 
effects are quite small since 4 He is hardly excited. We 
can then expect that the scattering is described by the 
optical potential that is obtained by folding the local ver¬ 
sion of microscopic £/,;x with the 4 He density. We refer 
to this model as the double single-folding (DSF) model 
in this paper. Very recently, it was shown that the DSF 
model well accounts for differential elastic cross sections 
and ctr for 4 He scattering without introducing any ad¬ 
justable parameter [36]. In fact, the DSF model yields 
better agreement with the data than the DF-FDA model 
for 4 He scattering. If the success of the model Hamilto¬ 
nian (2) for deuteron and 4 He scattering is not accidental, 
the model Hamiltonian should be good also for 3 He scat¬ 
tering. This is an important question to understand AA 
scattering systematically. 

In this paper, we investigate how good the model 
Hamiltonian is for 3 He scattering, and show that the 
model Hamiltonian works well and the DSF model based 
on the model Hamiltonian is a practical model to de¬ 
scribe 3 He scattering. This analysis is made in a wide 
incident-energy range of 30-150 MeV/nucleon. We con¬ 
sider heavier targets such as 58 Ni and 208 Pb in order to 
make our discussion clear, since the ^-matrix is evaluated 
in nuclear matter and hence the ^-matrix folding model is 
considered to be more reliable for heavier targets. 3 He is 
more fragile than 4 He, but less fragile than deuteron. We 
then investigate three nucleon breakup of 3 He by using 
CDCC. In the model Hamiltonian (2), the nucleon-target 


potential t/,;x is local. This makes CDCC calculations 
feasible. We show that the effects are appreciable only 
at lower incident energies around 40 MeV/nucleon and 
negligibly small at higher incident energies. This makes 
the DSF model reliable. In fact, the model yields better 
agreement with the experimental data on total reaction 
and differential cross sections than the DF-FDA model. 
We also investigate how the spin-orbit force between 3 He 
and T affects 3 He scattering by using the DSF model. Fi¬ 
nally we analyze effects of knock-on exchange processes 
on 3 He scattering. For AA scattering, the processes make 
the microscopic optical potential nonlocal, but the pro¬ 
cesses are approximately treated in the DSF model, since 
Ui t is localized. We then investigate how the approxi¬ 
mation affects 3 He scattering. 

In Sec. II, we derive the model Hamiltonian (2) from 
the many-body Hamiltonian with reasonable approxima¬ 
tions, using the multiple scattering theory [27-29]. Brief 
explanation is made on the DSF and DF-FDA models. 
The explicit form of the spin-orbit potential between 3 He 
and T is shown, and four-body CDCC is recapitulated. 
In Sec. Ill, numerical results are shown. Section IV is 
devoted to a summary. 


II. MODEL BUILDING 


We consider the scattering of P with mass number Ap 
from T with mass number Ax- In principle, the scatter¬ 
ing is described by the many-body Schrodinger equation 


Kr + hp + h'p + 


v v - E 

iePpeT 


T (+) = 0 


(3) 


for the total wave function where Vij is the re¬ 

alistic NN interaction and hr stands for the internal 
Hamiltonian of T. The total energy E is related to the 
incident energy E H 11 in the center of mass system as 
E = Ef™ + e 0 (P) + eo(T), where e 0 (P) and eo(T) are 
the ground-state energies of P and T, respectively. Fol¬ 
lowing the multiple scattering theory [27-29], one can 
rewrite Eq. (3) into 


Kr + hp + /lx + 


Y - 1 
Y 


J2 r b 

ieP.jeT 


E 




0, (4) 


where t 13 denotes the effective NN interaction in nuclear 
medium and the factor Y = ApA^p represents the number 
of the Tij working between P and T. For A A scattering 
with Y 1, the factor (Y — 1 )/Y can be approximated 
into 1. When Eq. (4) is derived from Eq. (3), the an- 
tisymmetrization between nucleons in P and those in T 
are neglected. However, the antisynnnetrization effects 
are well taken care of, if the T,y are symmetrical with re¬ 
spect to the exchange of colliding nucleons [37, 38]. The 
is often replaced by the g -matrix (gij) in many appli¬ 
cations, since both include nuclear medium effects. 
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Therefore, we reach the Schrodinger equation 

K r + h P + h T + 9ij - E \ ^ (+) = 0. (5) 

ieP.jeT 

In this work, we use the Melbourne g-matrix [11] as gij. 
As mentioned in Sec. I, g(p) is evaluated in nuclear mat¬ 
ter by solving nucleon scattering on a single Fermi sphere. 
For consistency with the nuclear-matter calculation, we 
consider the nucleon-target subsystem in the P+T sys¬ 
tem and assume 

P = pT{r m ) ( 6 ) 

as p in g(p). This procedure is referred to as the target- 
density approximation (TDA) in this paper. 

In the TDA, g(pp) includes target-excitation effects 
approximately, but does not include projectile-excitation 
effects. We can then assume that »]>(+) = </>o(T)^, where 
(j>o (T) is the ground state of T and 'ip describes the scat¬ 
tering of P from T in its ground state. Left-multiplying 
Eq. (5) by </>o(T), we can get the Schrodinger equation 

[H eS - E?™ - e 0 (P)] = 0 (7) 

for ip. Here the nucleon-target potential 

t/ iT HMT)|£fe(p T )|<A>(T)> (8) 

i€T 

is composed of the direct and knock-on exchange terms. 
The knock-on exchange process makes Ui t nonlocal, but 
it can be localized with the Brieva-Rook method [6] 
based on the local semi-classical approximation with 
good accuracy [20]. In this approach, projectile excita¬ 
tions have to be treated explicitly by solving Eq. (7). 

When projectile excitations are negligible, one can as¬ 
sume ^p = (po{P)xo{R) with the ground state </>o(P) of 
P and the relative wave function xo{R)- Equation (7) is 
then reduced to the Schrodinger equation 

[Kr + Up>sf(R) - E™]xo(R) = o, (9) 

where the optical potential Hdsf is obtained by 


then obtained from Eq. (5) as 

tC A (#) = (<I>o| £ 9i j{pp + Pt)|3>o) (11) 

ieP.jeT 

for $0 = 0o(P)</>o(T). The potential is also composed of 
the direct and knock-on exchange terms. The nonlocal¬ 
ity coming from the knock-on exchange process can be 
localized by the Brieva-Rook method with good accu¬ 
racy [21]. In actual calculations, /7c 0 ui is added to the 
localized U^,p A . This is the DF-FDA model. 

Comparing Eq. (10) with Eq. (11), one can see that 
the difference between the DSF and DF-FDA potentials 
mainly comes from that between the TDA and the FDA; 
see Appendix A for the detail. Since pp + pp > Pt , 
nuclear-medium effects are smaller in the DSF potential 
than in the DF-FDA potential. This makes the DSF 
potential more attractive and more absorptive than the 
DF-FDA potential. 

As for the central part of DSF and DF potentials, we 
summarize the explicit forms in Appendix A. In general, 
the 3 He optical potential has the spin-orbit part in addi¬ 
tion to the central part. General derivation of the spin- 
orbit part for AA scattering was shown in Ref. [39]. For 
3 He scattering, the spin-orbit part was calculated in Ref. 
[40, 41] for the DF model and in Ref. [41] for the SF 
model. Although adjustable parameters are introduced 
in the analyses, spin-orbit effects on differential cross sec¬ 
tions are similar between the two models. We then derive 
the spin-orbit part of 3 He optical potential only in the 
DSF model. In 3 He, the two-proton subsystem is con¬ 
sidered to be spin-singlet with good accuracy. Hence the 
spin-orbit force U^ F (R)L ■ s of 3 He scattering at an in¬ 
cident energy E- ln in the laboratory system is obtained by 
folding the spin-orbit force U^(r)i ■ s of neutron scat¬ 
tering at an incident energy £j n /3 with the neutron wave 
function <p n (t) in 3 He, where s is the neutron spin, r ( t ) 
is the coordinate of neutron from the center of mass of T 
( 3 He) and L (£) is the angular momentum with respect 
to R (r). Assuming that (p n (t) is an s-state, one can get 

Uv%p(R)L-s = J cp n (t)*U A p(r)l ■ scp n (t)dt, 


Uu SF (R) = (MP)\ U ™\MP))- ( 10 ) 

ieP 


Ap + Ap 1 dZ{R) L 
A P (1 + Ap) R dR 


( 12 ) 


The potential 1/dsf is local, since Utp is localized. This 
folding procedure is the DSF. The difference between the 
DSF and SF models comes from the difference of Uip. In 
the DSF model, Uip is a microscopic potential obtained 
by folding the Melbourne g- matrix with pp, but Utp is 
a phenomenological nucleon-nucleus potential in the SF 
model. The Coulomb potential t/coui is added to Hdsf 
in actual calculations. It is reported in Ref. [36] that the 
DSF model works well for 4 He scattering. 

Another approach to AA scattering is the DF-FDA 
model. In this approach, the g-matrix g{pp + pp) in¬ 
cludes both projectile- and target-excitation effects ap¬ 
proximately. The optical potential of P+T scattering is 


with 

Z(R) = J e iK R ^^p-p(K)dK, (13) 

U(K) = J e~ iK r U^(r)dr, (14) 

P{K) = J e lK t \(p n {t)\ 2 dt, (15) 

where the explicit form of U^p{r) is shown in Ref. [12]. 

Finally, three-body breakup effects of 3 He are taken 
into account with four-body CDCC [3]. Four-body dy¬ 
namics of the N+N+N +T system is treated in a model 
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space 

/ ~lm 

V = Y, |0 7 (P))(^(P)|, (16) 

7=0 

where </ 7 (P) is the 7 th eigenstate with an eigenenergy 
e 7 (P) obtained by diagonalizing hp by the Gaussian basis 
functions and among the eigenstates the 7 m th eigenstate 
<t > 7m (P) has the highest eigenenergy in the model space V; 
note that P is 3 He in the present case. This model-space 
approximation reduces Eq. (7) to 

V[H eS - E™ - e 0 (P m = 0, (17) 

where 

7 m 

P'ip = £^r(P )x,(R). (18) 

7=0 

This leads to the CDCC equation for x 7 (fi) as 
[E i c n m -Afl-(e 7 (P)- eo (P))] X 7 ( J R) 

7 m 

= 5> 7 (P)lE U i T |<MP))X7 '(«)• (19) 

7 ' ieP 

Here the folding potential U,n and the coupling potentials 
(^ 7 (P)| X)iep UiT |^ 7 / (P)) are calculated with direct nu¬ 
merical integration in the coordinate space. The CDCC 
equation is numerically solved with kinematical relativis¬ 
tic corrections, but the corrections are negligible in the 
present cases. 

As shown later in Sec. Ill, 3 He breakup effects are 
appreciable only at lower incident energies around 40 
MeV/nucleon and negligible at higher incident energies. 
Therefore, we estimate the breakup effects as simply as 
possible. For this purpose, we assume that 3 He is a 
spinless particle, and take the Minnesota force [42] as 
the nucleon-nucleon interaction in 3 He and introduce a 
three-body force to reproduce the binding energy of 3 He. 
In CDCC calculations, 0 + and 2 + states are considered 
as breakup states, since the contribution of 1 “ breakup 
states is confirmed to be negligibly small. For each of 0 + 
and 2 + , we take 35 breakup states below 20 MeV. The 
model space spanned by the ground and breakup states 
is confirmed to give the convergence of CDCC solutions 
for the elastic and reaction cross sections. 


III. RESULTS 

We consider 3 He scattering from 58 Ni and 208 Pb tar¬ 
gets in a wide incident-energy range of 30 < E- ln /Ap < 
150 MeV, where E- In is the incident energy in the lab¬ 
oratory system. First the scattering are analyzed with 
the DSF and DF-FDA models in which the spin-orbit 
force is not included. As the 3 He density (pp), we con¬ 
sider the density calculated with the three-nucleon model 


mentioned in Sec. II and the phenomenological density 
determined from electron scattering [43]. For the latter, 
finite-size effects due to the proton charge are unfolded 
in the standard manner [44], and the neutron density is 
assumed to have the same geometry as the proton one. 
Both the densities yield almost the same differential cross 
section, so we will take the former density to compare re¬ 
sults of the DSF model with those of CDCC later. 

The target density px is calculated with the spherical 
Hartree-Fock (HF) method with the Gogny-DlS interac¬ 
tion [45] in which the spurious center-of-mass motion is 
removed with the standard procedure [19]. In the present 
calculation, pairing effects are not included, but it is 
possible to take into account the effects with the spheri¬ 
cal Hartree-Fock-Bogoliubov (HFB) method with Gogny- 
DlS force. We have confirmed that the effects are negli¬ 
gible on ctr and differential cross sections at q < 5 fm _1 
for 3 He+ 58 Ni scattering at E- In /Ap = 72 MeV, where q 
is the transfer momentum. We can also deal with effects 
of nuclear deformation and inelastic channels by solv¬ 
ing the coupled-channel equation with coupling poten¬ 
tials among elastic and inelastic channels. The coupling 
potentials can be obtained by folding g-matrix with the 
transition densities calculated by structure models; see 
for example Refs. [19, 46, 47]. 

The Melbourne interaction is provided only up to kp = 
1.5 fuU 1 (p = 1.37po), where kp is the Fermi momentum 
and po is the normal density. We then assume that the 
Melbourne interaction at kp > 1.5 fm _1 is the same as 
that at kp = 1.5 fm _1 . This assumption does not affect 
any result of DSF calculations, since the g-matrix at kp > 
1.5 fm _1 yields small effects only on differential cross 
sections at larger angles in DF-FDA calculations. 

Figure 1 shows ctr, for 3 He scattering from 58 Ni and 
208 Pb targets in a range of E[ n /Ap = 30-150 MeV. 
The DSF model (circles) yields better agreement with 
the experimental data [48, 49] than the DF-FDA model 
(squares). For 58 Ni target, the DSF model slightly under¬ 
estimates the data around E- ln /Ap = 30 MeV, but this 
underestimation is solved by projectile breakup effects, 
as shown later in Fig. 3(b). 

Precisely speaking, the results of Fig. 1 depend on 
which g-matrix is taken. However, it is shown in Ref. 
[ 11 ] that elastic scattering are mainly determined by the 
on-shell part of g-matrix. The Melbourne g-matrix is 
constructed from Bonn-B nuclear force [50]. We have 
confirmed that the g-matrix constructed from CD-Bonn 
nuclear force [51] yields almost the same results as the 
Melbourne g-matrix; for example, the difference is only 5 
mb for 3 He+ 208 Pb scattering at E- m /Ap = 72.3 MeV. We 
will make further discuss on this point in the forthcoming 
paper. 

Figure 2 shows differential cross sections da/dtt as a 
function of transfer momentum q for 58 Ni and 208 Pb tar¬ 
gets. For Ei n /Ap ss 40 and 150 MeV, the DSF model 
(solid line) definitely yields better agreement with the ex¬ 
perimental data [52-56] than the DF-FDA model (dashed 
line). For E- In /Ap = 72 MeV, agreement with the data 
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FIG. 1: (Color online) Total reaction cross section <tr as 
a function of E\ n /Ap for (a) 3 He+ 58 Ni scattering and (b) 
3 He+ 208 Pb scattering. The circles (squares) stand for results 
of the DSF (DF-FDA) model. The spin-orbit force is not in¬ 
cluded in both the models. The experimental data are taken 
from [48, 49]. 


is comparable between the two models; more precisely, 
the DSF model is better than the DF-FDA model at 
q < 2 fm - , but the latter is superior to the former at 
q > 2 fm -1 . The overestimation of the DSF model at 
q > 2 fm -1 comes from the fact that knock-on processes 
are treated only approximately in the DSF model, as dis¬ 
cussed later in Fig. 5. The DSF model thus yield better 
description for 3 He scattering than the DF-FDA model. 

Figure 3 shows projectile-breakup and spin-orbit force 
effects on da/dtt and or for 3 He+ 58 Ni scattering. In 
panel (a) for da/dfl, the dot-dashed and dashed lines 
stand for results of the DSF model with and without 
the spin-orbit force respectively, while the solid lines de¬ 
note results of CDCC in which the spin-orbit force is ne¬ 
glected. Hence the difference between the dashed and 
solid lines shows projectile-breakup effects, while that 
between the dashed and dot-dashed lines corresponds 
to spin-orbit force effects. Both the effects are small, 
but improve agreement with the data. More precisely, 
spin-orbit force effects are appreciable at higher inci¬ 
dent energies around E- m /Ap = 150 MeV, but projectile- 
breakup effects are visible at lower incident energies such 
as E' ln jAp = 40 MeV. The present result on spin-orbit 
force effects is consistent with the previous one of Ref. 
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FIG. 2: (Color online) Differential cross sections as a function 
of transfer momentum q for (a) 3 He+ 58 Ni scattering and (b) 
3 He+ 208 Pb scattering. The cross section at each E ln /Ap is 
multiplied by the factor shown in the panel. The solid and 
dashed lines denote results of the DSF and DF-FDA models, 
respectively. The spin-orbit force is not included in both the 
models. The experimental data are taken from Refs. [52-56]. 


[41]. 

In panel (b) for or, the circles denote results of the 
DSF model without the spin-orbit force, while the trian¬ 
gles correspond to results of CDCC in which the spin- 
orbit force is neglected. In this panel, results of the DSF 
model with the spin-orbit force are not shown, since they 
agree with those of the DSF model without the spin- 
orbit force. This indicates that spin-orbit force effects are 



















6 


negligibly small. Projectile-breakup effects are shown by 
the difference between the triangle and the corresponding 
circle. The effects are more appreciable for ctr than for 
dcr/dfl. The effects improve agreement with the experi¬ 
mental data particularly at low incident energies around 
Ein/Ap = 30 MeV. 

Comparing Fig. 3 with Figs. 1 and 2, one can easily 
see that CDCC calculations based on the model Hamilto¬ 
nian well describe 3 He scattering and yield better agree¬ 
ment with the data than the DF-FDA model. The model 
Hamiltonian is thus good not only for deuteron and 4 He 
scattering but also for 3 He scattering. 

Figure 4 shows L dependence of the elastic S'-matrix 
elements Sl- The filled (open) circles denote the Sl cal¬ 
culated with the DSF model (CDCC), where the spin- 
orbit force is neglected. Projectile-breakup effects shown 
by the difference between the open and filled circles are 
sizable for large L but not for small L. At large L , fur¬ 
thermore, the effects reduce the absolute value of Sl. 
Therefore the dynamical polarization potential generated 
by 3 He breakup is strongly absorptive in the peripheral 
region of T. This is the reason why 3 He breakup effects 
are more appreciable for or than for da/dfl. 

In the DSF model, the nucleon-target potential Ui t is 
localized by the Brieva-Rook method. Now we discuss 
how the localization affects 3 He scattering, although the 
localization is accurate for nucleon-target scattering itself 
as shown in Ref. [20]. When the nucleon-target poten¬ 
tial is not localized, the potential between P and T is 
obtained by 

^ F DA = <*o| E S«G*r)l*o>. (20) 

iePpeT 


This is nothing but the DF model with the TDA. This 
model is referred to as the DF-TDA model in this paper. 
The explicit form of DF-TDA potential is shown in Ap¬ 
pendix A, together with the explicit forms of DF-FDA 
and DSF potentials. 

In the DF-TDA model, the nonlocality of the folding 
potential comes from knock-on exchange processes be¬ 
tween interacting two nucleons, one in P and the other in 
T. As a consequence of the exchange, the one-body densi¬ 
ties pp(rp) and pt(j’t) are changed into the correspond¬ 
ing mixed densities /5p(rp, rp — s) and /5 t(^t, t’t + s ) as 
shown in Appendix A, where rp (rx) denotes the coor¬ 
dinate of nucleon in P (T) from the center of mass of P 
(T) and s is the relative coordinate between interacting 
two nucleons. The one-body and mixed densities of P 
are defined by 

PpOp) = E^pE rp ) < Ea( rp )’ ( 21 ) 

a 

Pp(rp,r F - s) = E ^E( rp )^P;a( rp “ s ) ( 22 ) 
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FIG. 3: (Color online) Projectile-breakup and spin-orbit force 
effects on (a) differential cross sections and (b) total reaction 
cross sections for 3 He+ 58 Ni scattering. In panel (a), the dif¬ 
ferential cross section at each Ei n /Ap is multiplied by the 
factor shown in the panel. The dot-dashed and dashed lines 
denote results of the DSF model with and without the spin- 
orbit force, respectively, while the solid lines correspond to 
results of CDCC in which the spin-orbit force is neglected. 
In panel (b), the circles stand for results of the DSF model 
without the spin-orbit force, while the triangles correspond 
to results of CDCC in which the spin-orbit force is neglected. 
Results of the DSF model with the spin-orbit force agree with 
those of the DSF model without the spin-orbit force. The ex¬ 
perimental data are taken from Refs. [48, 49, 52-54], For 
Ei n /Ap = 72.3 MeV in panel (a), the CDCC result (solid 
line) is somewhat different from the DSF one (dashed line) at 
1.5 < q < 2.5 fin -1 , although the two results are close to each 
other at backward angles q > 2.5 fin -1 ; note that the dashed 
line agrees with the dot-dashed line at 1.5 < q < 2.5 fm -1 . 
A few percent deviation between CDCC and DSF results in 
panel (b) for or mainly comes from the difference between 
CDCC and DSF results in panel (a) at 1.5 < q < 2.5 fm -1 . 


with the single particle wave function 4>p. a {rp) of P char¬ 
acterized with the quantum number a, where fj, denotes 










7 



FIG. 4: (Color online) Projectile-breakup effects on elas¬ 
tic S-matrix elements for 3 He+ 58 Ni scattering at Ei u /Ap = 
40 MeV. The filled (open) circles correspond the elastic S- 
matrix elements calculated with the DSF model (CDCC); 
here note that the spin-orbit force is not included. The elastic 
S-matrix elements tend to a point (1,0) as L increases. 


either proton or neutron. One can take the same defini¬ 
tion also for the one-body and mixed densities of T. 

The knock-on exchange is properly treated in the DF- 
TDA model. The resulting nonlocal potential between P 
and T can be localized with high accuracy by the Brieva- 
Rook method based on the local semi-classical approxi¬ 
mation [21]. In the DSF model, the change of pt (f t) 
due to knock-on exchange is taken into account in the 
nucleon-target potential. However, the change of pp(rp) 
is not treated as shown in Appendix A, since the nucleon- 
target potential is localized. Therefore, one can see how 
the change of pp(r’p) due to knock-on exchange affects 
3 He scattering, comparing results of the DSF model with 
those of the DF-TDA model. 

Figure 5 shows differential cross sections calculated 
with the DSF, DF-TDA and DF-FDA models for (a) 
3 He+ 58 Ni scattering and (b) 3 He+ 208 Pb scattering, 
where the spin-orbit force is neglected. The dashed 
(solid) lines correspond to results of the DF-TDA (DSF) 
model. The difference between the dashed and solid lines 
shows effects of the pp change due to knock-on exchange. 
The effects are small for higher incident energies such as 
E-m/Ap = 150 MeV. For lower incident energies less than 
E-m/Ap = 100 MeV, meanwhile, the effects appear at 
larger angles and improve agreement with the experimen¬ 
tal data, particularly in q > 3 fm -1 for 3 He+ 58 Ni scatter¬ 
ing at E[ n /Ap = 40 MeV and in q > 2 fm -1 for 3 He+ 58 Ni 
scattering at E- m /Ap = 72 MeV. The DF-TDA model 
thus yields better agreement with the data than the DSF 
model particularly at the larger angles. Eventually, the 
DF-TDA model (dashed line) well reproduce the data 
in a wide incident-energy range from E- m /Ap = 40 MeV 
to 150 MeV, compared with the DF-FDA model (dot- 
dashed line). For <tr, furthermore, we have confirmed 
that the DF-TDA model yields almost the same result 
as the DSF model and hence better agreement with the 


data than DF-FDA model. 



q [fm' 1 ] 



q [fm' 1 ] 


FIG. 5: (Color online) Comparison of DSF, DF-TDA and DF- 
FDA models for differential cross sections in (a) 3 He+ 58 Ni 
scattering and (b) 3 He+ 208 Pb scattering. The differential 
cross section at each -EA/Ap is multiplied by the factor shown 
in the panel. The dashed (dot-dashed) lines denote results of 
the DF-TDA (DF-FDA) model, while the solid lines stand for 
results of the DSF model. The spin-orbit force is not included 
in both the models. The experimental data are taken from 
Refs. [52-56]. 

We compare DSF, DF-TDA and DF-FDA potentials 
in Fig. 6 for 3 He+ 58 Ni scattering at E- m /Ap = 40 MeV. 
The FDA has stronger Pauli-blocking effects than the 
TDA because of pp + pt > Pt, so that the DF-FDA 
potential (dot-dashed line) is less attractive and less ab¬ 
sorptive than the DF-TDA potential (dashed line). The 

















DSF model (solid line) well simulates the DF-TDA po¬ 
tential in the peripheral region. This is the reason why 
the DSF model well simulates the DF-TDA model at for¬ 
ward angles of q < 3 fm -1 . 

Finally, we comment on a traditional Woods-Saxon 
(WS) parametrization for the phenomenological optical 
potential. For 3 He scattering, it is difficult to reproduce 
the folding potential with a single WS form in the full 
range of R. In general, the forward scattering is sensi¬ 
tive to the peripheral region of the folding potential. If a 
single WS form is determined from the folding potential 
in the peripheral region, the single WS potential can re¬ 
produce the differential cross section at q < 2 fm -1 . If a 
double WS form is used, the potential is close to the fold¬ 
ing potential and reproduces the differential cross section 
up to q = 4 fm -1 . Thus the double WS form is recom- 
mendable for the present scattering. 




FIG. 6: (Color online) Microscopic optical potentials for 
3 He+ 58 Ni scattering at Ein/Ap = 40 MeV. Here V(R) 
(W(R)) is the real (imaginary) part of the potentials. The 
dashed (dot-dashed) line stands for the DF-TDA (DF-FDA) 
potential, whereas the solid line denotes the DSF potential. 


IV. SUMMARY 

We have investigated how good the model Hamilto¬ 
nian (2) is for 3 He scattering. In the model Hamiltonian, 
the nucleon-target potential [/,t is obtained by folding 
the Melbourne g-matrix with the target density and lo¬ 
calizing the resultant nonlocal folding potential with the 
Brieva- Rook method. In many calculations done so far, 
the phenomenological nucleon optical potential was used 
as Ui t- As an advantage of the present approach from 
the previous one, the microscopic Ui t is obtainable even 
for the nucleon+target system in which no experimental 
data is available. As an advantage of the present ap¬ 
proach from the DF-FDA model, the model Hamiltonian 
is consistent with the fact that the single Fermi sphere 
is considered when the ^-matrix is evaluated in nuclear 
matter. In the present approach, projectile excitation ef¬ 
fects have to be treated explicitly, but this can be done 
with CDCC since UiT is localized. 

The validity of the model Hamiltonian has been in¬ 
vestigated in a wide incident-energy range of Ej n /Ap = 
30-150 MeV for heavier targets 58 Ni and 208 Pb, since the 
^-matrix is considered to be more reliable for heavier tar¬ 
gets. CDCC calculations well reproduce the experimen¬ 
tal data with no adjustable parameter. The calculations 
show that 3 He breakup effects are appreciable only at 
lower incident energies around 40 MeV/nucleon and neg¬ 
ligibly small at higher incident energies. Therefore the 
DSF model also well accounts for the experimental data. 
In fact, the DSF model yields better agreement with the 
data than the DF-FDA model. Thus we can propose that 
the DSF model is a practical model to describe 3 He scat¬ 
tering and can conclude that the model Hamiltonian is 
good not only for deuteron and 4 He scattering but also 
for 3 He scattering. 

As a merit of the DSF model, the spin-orbit force be¬ 
tween P and T can be easily evaluated. Effects of the 
force are found to be appreciable only at higher incident 
energies such as Ein/Ap = 150 MeV and negligibly small 
at lower incident energies. This is consistent with the 
previous analysis of Ref. [41] based on the phenomeno¬ 
logical nucleon-target potential. 

When the nucleon-target potential is not localized, 
the DSF model becomes the DF-TDA model. Knock- 
on exchange processes are properly taken into account 
in the DF-TDA model, but the change of pp(rp) due 
to knock-on exchange is not in the DSF model because 
the nucleon-target potential is localized. For lower in¬ 
cident energies such as Ej n /Ap < 70 MeV, effects of 
the pp change appear in differential cross sections at 
larger angles of q > 2 fnV 1 and improves agreement 
with the experimental data there. When one is inter¬ 
ested in the backward angles, the DF-TDA model is bet¬ 
ter than DSF model. For higher incident energies such as 
Ei n /Ap = 150 MeV, meanwhile, the effects are small and 
hence the DSF model is as good as the DF-TDA model. 

In conclusion, the model Hamiltonian (2) works well 
for scattering of Os-shell nuclei such as deuteron, 3 He 
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and 4 He. As a future work, it is quite interesting to 
investigate whether the model Hamiltonian is good also 
for scattering of Op-shell nuclei such as 12 C and 16 O. 
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Appendix A: Explicit forms of DF-FDA, DF-TDA 
and DSF potentials 

Here we show the explicit forms of DF-FDA, DF-TDA 
and DSF potentials and discuss their relation. This is an 
overview of the previous discussion in Ref. [36]. 

The DF potential Hdf for P+T scattering at an inci¬ 
dent energy E- m is composed of the direct and knock-on 
exchange terms, U^p and U% p, defined by [15, 18, 19] 

USf(R) = Y1 [ Pp( r p)pT(rT)g™(s;p^)drpdr T ,{Al) 

UBf(R) = Y1 [ Pp{rp,rp - s)p^{r T} r T + s) 

ll,V ^ 

x 3^0; P^) ex P [-iK(R) ■ s/M]drpdr T , (A2) 

where s = rp — rx + R and p and v denote either pro¬ 
ton or neutron. The non-local U EX has been localized 
with the local semi-classical approximation [6], where the 
local momentum HK(R) between P and T is given by 
HK(R ) = yj2M(E- m — Uuf(R)) with the reduced mass 
M = ApA-p/(Ap + At) of the P+T system. The direct 
and exchange parts, g EE and g EX , of g-matrix depend on 
the density p defined by 

Pnv = Pp( r P — s/2) + Px( r t + s /2) (A3) 

in the FDA and 

Pfj.u = Pt(t t + s/2) (A4) 

in the TDA. The difference between the DF-FDA and 
DF-TDA stems from only that between Eqs. (A3) and 
(A4). 

For N+T scattering at an incident energy E^, the SF 
potentials for proton (p = —1/2) and neutron (p = 
1 /2) scattering also consist of the direct and knock-on 


exchange terms, U ER and U EX , defined by [20] 


u B R (. r v) = Yl [ PT(rT)g™(s;Puv)dr T , (A5) 

V J 

u » X ( r n) = '52 [ PT(rT,r T + s) 

X /v) exp l-iK^r^) • s]cir T (A 6 ) 


for s = r M — t*T) where r M is the coordinate of an in¬ 
cident nucleon N from the center of mass of T and the 
local momentum hK ^(r M ) between N and T is given by 
hK^ir^) = sj2pw'r(Ef n - C/ M (r M )) for the reduced mass 
Pnt of the N+T system. 

The DSF potential C/dsf is obtained by folding the 
nucleon-target potentials U ER + U EX with the projectile 
density p 1 /,. The direct and knock-on exchange parts are 
then given by 


1®(H)=£ / p£(r P )UB R (R + rp)drp 

= £ / Pp( r p)PT(rT)g^ R {s; p^)drpdr T , (A7) 

(L,l/ 

UBsf(R) = Y1 [ Pp{rp)Uf^{R + r P )dr P 

= £ / Pp( r p)PT( r T,r T + s) 

X Pnv) ex P [~ iK n( r n) • s/M]drpdr T , (A 8 ) 


Now we consider heavy targets satisfying At Ap > 
1 for simplicity. The g matrix depends on an energy 
of nucleon in P. The energy is E x t for N+T scattering 
and E- ln /Ap for P+T scattering. We can then find that 
^DF-TDA = ^DSF when E hi = E in/Ap■ In the periph¬ 
eral region of T important for elastic scattering, the local 
momenta /(R' M (r /i ) and HK(R) are close to their asymp¬ 
totic values, hKfj,( oo) and HK(oo), respectively. When 
E = Ei n /Ap, the asymptotic values satisfy 


K lt (oo) = K(oo)/M. (A9) 


Therefore we see that U// 'f-tda = ^dsf if Pp(rp, rp — s) 
is identical with /Op(rp). The approximation pp(rp, rp — 
s) w Pp(rp,rp) = pp(rp) is known to be good in the 
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peripheral region of T [20]. 
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